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GLOSSARY 


Hg 

Mercury 

ICAO 

International  Civil  Aviation  Organization 

IFF 

Identification  Friend  or  Foe 

kft 

Kilofeet 

mb 

Millibars 

MESO 

Military  Electronic  Systems  Operations 

MSE 

Mean -Square-Error 

N 

North 
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National  Bureau  of  Standards 
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Nautical  Mile(s) 
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Root-Mean  -Square 

SECTIO'  T 


INTRODUCTION 


Buyers  of  modern  Air  Defense  radars  usually  specify  bounds  for  the  total  root-moan- 
square  (rms)  errors  of  range,  azimuth  and  height.  These  errors  have  bias  and  random 
(jitter  and  thermal)  components  which  appear  as  summed  variances.  With  a  colocated  pre¬ 
cision  track  radar  with  known  errors  much  smaller  than  the  radar  under  test,  the  total  rms 
error  Is  estimated  in  a  straightforward  manner.  *  For  test  sites  without  such  a  standard, 
the  problem  is  more  difficult.  Such  is  the  case  for  the  majority  of  buyers  of  the  GE-592 
3-D  Air  Defense  Radar.  However,  a  standardized  approaoh  to  this  problem  has  emerged 
at  Military  Electronic  Systems  Operations  (MESO)  in  which  the  bias  and  random  errors  are 
characterized  separately,  without  the  need  for  a  colocated  precision  track  radar.  A  methc  1 
of  making  unbiased  estimates  of  the  variance  of  random  radar  errors  has  been  developed 
by  this  author  and  a  paper  is  currently  In  preparation.  This  paper  deals  entirely  with  the 
characterization  of  radar  height  bias  errors. 


*  The  requirement  for  colocation  stems  from  a  need  for  accurate  time  registration  between 
the  two  sites  and  to  eliminate  disparate  environmental  effeits  from  olutter,  multipath, 
and  propagation. 


1-1/1-2 


SECTION  n 


STATISTICAL  DESCRIPTION  OF  RADAR  ERRORS 


For  all  practical  purposes,  radar  errors  can  be  classified  as  either  fixed  or  random. 
The  fixed  errors  typically  result  from  an  error  in  design  implementation  (hardware  or  soft¬ 
ware)  for  which  beam  misplacement  or  pedestal  tilt  are  examples.  Random  errors  have 
various  types  of  probability  densities  (e.g,  uniform  and  Gaussian)  and  can  be  zero-mean  or 
have  a  finite  mean.  A  nonzero  mean  can  result  from  a  calibration  error  for  instance,  or 
simply  be  too  costly  to  entirely  eliminate.  Often  this  mean  may  vary  systematically  with 
time.  Uniform  probability  densities  arise  from  sources  like  quantization.  Gaussian  den¬ 
sities  arise  from  a  collection  of  sources  like  equipment-induced  Jitter  for  example.  From 
a  sampling  standpoint,  an  important  consideration  is  the  correlation  time  from  sample  to 
sample.  Many  errors  are  independent  between  samples.  Such  is  the  case  for  thermal  and 
jitter  errors  which  have  a  correlation  time  which  is  short  compared  to  a  revisit  period. 

Often  these  errors  are  simply  referred  to  as  "random".  But  many  bias  errors  are  random 
as  well,  in  that  samples  can  be  quite  different  and  unpredictable  from  time  to  time  although 
the  time  scale  for  which  they  seem  constant  may  be  quite  large;  hours  or  even  days.  Solar, 
ambient  temperature,  and  average  environmental  factors  (like  refraction)  are  examples  of 
sources  of  variable  bias  errors  which  are  unpredictable  and  long  term  random  beyoud  the 
point  for  which  the  radar  designers  can  account  for  them  by  calibration  or  compensation. 
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ALTIMETER  ERRORS 


Height  bias  errors  can  derive  from  refraction,  electronic,  mechanical,  and  software 
sources.  In  the  absence  of  a  colocated  precision  track  radar,  a  convenient  source  of  a 
height  standard  is  the  identification  friend  or  foe  (IFF)  Mode-C  reply  from  the  interrogated 
target  aircraft.  This  IFF  height  comes  from  the  altimeter,  which  is  a  barometric  device 
employing  an  aviation  standard  international  civil  aviation  organization  (sometimes  called 
the  ICAO  standard)  atmosphere  model.  Unfortunately,  this  altimeter,  while  relatively 
pr  ,'C  -  r .  ^cin  be  highly  inaccurate  particularly  at  high  altitudes.  Some  discussion  of  this 
st a  <ava  uid  its  error  is  therefore  in  order. 


Figure  3-1  shows  the  ICAO  models  for  temperature  and  pressure.  Below  18  kft  the 
t'lot  scales  the  pressure  model  by  a  surface  value  obtained  from  the  nearest  weather 
taticn.  Abo.  j  18  kft  the  model  is  uncorrected  and  based  upon  a  surface  pressure  of  29.  92 
c.  lig  (mercury)  or  1013  mb  (millibars).  Often  the  pilot  will  use  an  ambient  temperature 
me  m  urement  to  further  compensate  his  altimeter  reading  with  the  aid  of  a  special  slide  rule. 
Fbwl  .3  is  ac- dally  doing  is  converting  the  pressure  altitude  into  a  corrected  density  altitude. 
First  the  barometer  in  the  altimeter  measures  the  pressure  p(H)  at  the  true  altitude  H. 

An  analog  computer  in  the  altimeter  then  enters  the  ICAO  standard  atmosphere  of  pressure 
versus  altitude  to  find  the  altimeter  setting  HQ(p).  This  may  be  an  analog  or  digital  process. 
Below  18  kft  this  setting  is  adjusted  upwards  approximately  1000  ft  for  every  1  in.  Hg 
(33. 857  mb)  that  the  nearest  applicable  measured  surface  pressure  exceeds  the  standard, 
i.  e.  29. 92  in.  Hg  (1013  mb).  The  height  value  shown  to  the  pilot  on  the  altimeter  is  thus 


Hc(kft) 


HQ(P)  + 


Pgjj  *  1013.  25 
36.922 


PgL  -  1013.  25 
269.  503 


(3-D 


for  barometric  pressures  less  than  500  mb  (IS  kft)  and  just  HQ(p)  for  all  altitudes  greater 
than  18  kft  (or  pressures  less  than  or  equal  to  500  mb).  The  level  transmitted  by  Mode-C 

4 1 

link  is  not  corrected  for  surface  pressure,  however,  and  is  just  HQ(p)  at  all  altitudes.  In 
Equation  (3-1),  the  following  definitions  apply: 

*  Except  perhaps  for  the  latest  military  aircraft  like  the  F-15. 


3-1 


H  -  Mode-C  height  (kft) 

c 

Hq  -  altimeter  reading  (kft) 

p  -  barometric  pressure  at  true  height  H,  (mb) 


PSL 


true  sea  level  pressure  (mb). 


There  are  three  principal  sources  of  uncertainty  in  the  Mode-C  transmitted  height. 

1.  Barometric  measurement  error 

2.  Mode-C  encoding  error 

3.  Natural  atmospheric  variability  from  the  standard  at  all  altitudes. 


3.1  MEASUREMENT  ERROR 

The  aircraft  altimeter  has  measurement  error  which  has  both  bias  and  random 
components.  There  doesn't  seem  to  be  any  data  available  on  the  bias  component,  but  the 
random  (zero-mean)  component  seems  well  characterized  by  a  standard  deviation  of  3  mb. 
The  surface  pressure  is  determined  by  smoothing  a  lot  of  data  and  has  a  measurement 
error  which  is  negligible. 


3.2  MODE-C  ENCODING  ERROR 

The  ICAO  standard  height  is  encoded  onto  the  Mode-C  transponder  link  with  a  95th 
percentile  error  of  125  ft.  This  implies  a  standard  deviation  of  0. 065  kft. 


3.3  ATMOSPHERIC  VARIABILITY 

Data  on  the  variability  of  pressure  at  sea  level  and  at  altitude  is  difficult  to  come  by. 
Two  approaches  are  followed  here.  First  consider  devlal  ons  from  standard  pressure  at 
various  north  latitudes  as  compiled  from  the  1966  U.  S.  Standard  Atmosphere  (l)e  and 
shown  in  Figure  3-2.  Presuming  a  uniform  distribution  between  these  limits  we  can  infer 
a  standard  deviation  and  normalize  by  the  mean.  The  results  are  plotted  in  Figure  3-3 
as  the  curve  connecting  circles  for  a  latitude  of  60°N.  Here  we  see  zero  variability  at  sea 
level,  an  obvious  contradiction  with  observation.  From  the  1976  U.S.  Standard  Atmosphere 
(2),  we  have  information  on  the  variability  of  temperature  and  density  shown  in  Figures  3-4 
and  3-5.  From  the  perfect  gas  law,  we  can  derive  an  expression  for  the  coefficient  of 
variation  pressure  in  terms  of  thn  coefficients  af  variation  for  temperature  and  density  and 
the  correlation  between  them.  Thus, 
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04.  ST  AND  AS  D  ATMOSPHERE  SUPPtEMiWTS, 


Figure  3-2.  Coefficient  of  Pressure  Variation  vs  Altitude 


PERCENT  DEPARTURE  PROM  1178  STANDARD  DENSITY 


Figure  3  -4.  Range  of  Systematic  Variability  of  Density  Around  the 
U.  S.  Standard  Atmosphere,  1976 


140  ISO  ISO  300  no  340  3S0  380  300  330 
TEMPERATURE, K 

Figure  3-5.  Range  of  Systematic  Variability  of  Temperature 
Around  the  U.  S.  Standard  Atmosphere,  1976. 
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<r  is  the  inferred  standard  deviatic  a  of  pressure  (mb), 

Q 

or  p  is  the  derived  standard  deviation  of  density  (gm/cm  ), 

<rT  is  the  derived  standard  deviation  temperature  (°K), 

p  _  is  the  correlation  coefficient  between  measured  density  and 
p  temperature,  i.e. 

ap  T^ap  °T  ’  and 

o 

p,  p,T  are  tha  pressure  (mb),  density  (gm/om  )  and  temperature  (°K) 
respectively. 

Because  we  don't  have  information  on  the  correlation  between  p  and  T  we  can  assume  a 
uniform  distribution  of  correlations  between  -1  and  1,  since  the  entire  range  1b  possible. 

In  other  words,  it's  just  as  likely  that  temperature  increases  while  density  decreases 
(increased  solar  heating  with  increasing  wind  velocity)  as  the  possibility  that  both  will  in¬ 
crease  simultaneously.  An  isothermal  process  is  also  possible  for  some  measurements, 
corresponding  to  p  pT=0  when  both  the  pressure  and  density  change  in  proportion. 

Integrating  Equation  (3-2)  over  all  possibilities  of  ppT  from  -1  to  1  eliminates  the  last 
te  m  which  is  equivalent  by  circumstance  to  zoro  correlation.  Thus, 

(3-3) 

From  Figures  3-4  and  3-5,  assuming  a  Gaussian  distribution  with  1%  tails  we  can  infer  stand¬ 
ard  deviations  for  density  and  temperature  as  input  for  Equation  (3-3).  The  results  are  plotted 
in  Figure  3-3  aB  the  triangles.  Since  the  direct  data  (circles)  about  a  / p  is  preferable 

Mk 

but  incorrect  at  low  altitudes,  we  use  the  inferred  data  at  low  altitudes  until  the  direct  data 
exceeds  it,  in  lieu  of  additional  information.  Thus  we  take  the  solid  line  of  Figure  3-3  as 
definitive  for  the  rest  of  the  repo:*t  and  interpret  it  ns  the  expected  variability  in  pressure 
from  standard  at  any  particular  altitude,  place,  time  of  day,  or  year  (presuming  ergodicity). 


i 

si 

\i 

[i 

[i  3. 4  ERROR  ANALYSE 

t  In  this  section  we  want  to  develop  a  reasonable  model  for  the  bias  error  by  using 

the  Mode-C  height  as  an  estimate  of  true  aircraft  altitude.  It  is  a  bias  error  in  the  sense 
that  the  erro .*  will  not  change  for  repeated  measurements  at  the  same  altitude  and  range  for 

S  times  short  compared  to  several  hours.  This  error  does  vary  as  a  function  of  target 

position  and  with  times  on  the  order  of  hours.  For  this  analysis,  we  assume  that  this  bias 
|  error  is  zero-mean  Gaussian  with  a  standard  deviation  that  can  be  estimated  from  the  data 
just  developed.  It  is  further  assumed  that  the  correlation  between  samples  is  on  the  order 
of  0. 9  for  samples  taken  15-minutes  apart  and  0. 1  for  samples  taken  4-hours  apart.  For 
simplicity  in  the  analysis  we  can  model  the  standard  atmosphere  as 

p(mb)  e  1033  exp  J  -  H-  ^— ■  J  (3-4) 

'  I  where 


a  =  204.  8  kft 

and 

b  =  (103.33  kft)2  . 

From  Equation  (3-1),  errors  in  the  barometric  altitude  and  encoder  induce  an  error  in  the 
Mode-C  height  of 

5H0  =  6H  4-  oHme  (3-6) 

where, 

<5H  *  the  transponder  encoding  error,  and  from  Equation  (3-4), 

vlQC 


I 

I 

j  I 


I 


2a  po 

<5  H  «=  - 

1  +  H/a 

which  relates  a  change  in  altimeter  pressure  to  a  change  in  standard  altitude  In  the  absenoe 
of  measurement  error.  Defining  6pm  as  the  measurement  error  and  <5p  ^  as  the  variation 
of  true  pressure  at  altitude  H. 
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Squaring  and  taking  expectation  values, 


o  2  2 

2  <7  +  <7 

Poo  P 


m 


H, 


(1  +  H/a)' 


+  cr  ■ 
H, 


enc 


where  from  par.  3. 2, 


(3-6) 


<r„  -  0.  065  kft. 

Henc 

The  standard  deviation  of  altimeter  error  is  plotted  in  Figure  3-6.  With  bias  errors  that 
have  standard  deviations  this  large,  the  altimeter  does  not  appear  to  be  a  reliable  standard 
of  reference  from  which  to  estimate  radar  height  bias  errors,  which  are  designed  to  be  far 
less  than  those  shown  for  the  altimeter.  Fortunately  we  can  correct  the  altimeter  after  the 
flight  by  the  careful  use  of  appropriate  radiosonde  data. 
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SECTION  IV 


I 


l 


i 


4 


RADIOSONDE 


4.1  DESCRIPTION  AND  HYDROSTATIC  CALCULATIONS 

Radiosonde  is  a  system  for  measuring  several  essential  properties  of  the  atmo¬ 
sphere;  static  pressure,  temperature,  relative  humidity,  and  wind  components  versus 
altitude.  The  instrumentation  is  carried  aloft  by  a  helium  balloon  and  the  measurements 
are  transmitted  by  radio  down  to  a  ground  station.  A  tracking  radar  is  used  to  estimate 
the  wind  components.  The  correct  altitude  is  estimated  on  the  ground  by  computer  from 
the  pressure,  temperature,  and  relative  humidity  measurements  by  integration  of  the 
hydrostatic  differential  equation 


(4-1) 


where 

p  is  ambient  pressure  (mb), 

H  is  altitude  (kft), 

a  -  10. 412243  "K/kft,  and 

Ty  Is  the  virtual  temperature  (#K), 

1.  e. , 

Ty  -  T(l-  0.379  pv)  (4-2) 

and 

T  is  the  ambient  dry-bulb  temperature  (*K), 

RH  is  the  relative  humidity,  and 

Py  is  the  vapor  pressure  at  temperature  T. 

The  vapor  pressure  is  the  partial  pressure  of  water  vapor  in  the  air  at  a  given  temperature  T 
It  is  an  empirical  function  given  in  Figure  4-1  along  with  reasonable  analytical  models. 

Note  that  the  virtual  temperature  reduces  to  the  ambient  temperature  when  the  relative 


humidity  is  zero. 
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.VAPOR  PRESSURE  OF  AIR  (mb) 


Two  common  approximations  made  for  the  temperature  as  a  function  of  altitude  are: 
(1)  constant  and  (2)  linear,  within  a  layer  defined  by  consecutive  radiosonde  transmissions.  * 
The  formulae  that  result  are  thus  useful  for  interpolating  between  radiosonde  reporting 
levels. 

For  a  constant-temperature  level,  ignoring  the  v-subscript  on  virtual  temperature, 


'■(*)•-■*  l 

H1 

H  -  Hi  +  T  <n  (t) 

P  -  Pj  <s’tPf‘%  (H-Hjlj 


where  the  subscript  1  refers  to  the  lower  altitude  report  for  the  layer  of  interest. 
For  a  variable  temperature  layer  (assumed  linear), 


Hi  V 


X  -  H, 


H-H*+  T>  [(t 


T2“T1 
a  (Hg-Hj) 


*  Sometimes  the  temperature  is  assumed  linear  with  pressure  as  a  means  of  defining 
the  layers.  In  this  case  the  hydrostatic  altitude  is, 

H  =  »l  -  IT  •  (VT0  +  [Ti-pi(^)]/n  (j)  • 
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or. 


-  Q 


1  + 


H  -  Hx 

v75? 


-  1 


H2-Hl 

VTX 


(4-6) 


where  the  subscript  2  refers  to  the  higher  altitude  report  for  the  layer  of  interest. 


4.2  POST-FLIGHT  ALTIMETER  CORRECTION  PROCEDURE 

The  post-flight  altimeter  correction  procedure  consists  of  the  following  steps: 

1.  As  a  preliminary  step,  append  virtual  temperature,  Equation  (4-2),  to 
the  radiosonde  data  table. 

2.  Infer  the  altimeter  barometric  pressure  from  the  ICAO  standard 
atmosphere  at  the  raw  altitude  computed  in  step  1. 

3.  Enter  the  radiosonde  data  table  in  the  layer  which  would  include  the 
pressure  inferred  in  step  2. 

4.  Compute  the  corrected  altitude  by  interpolation  according  to  Equations 
(4-3)  or  (4-5),  depending  upon  whether  the  virtual  temperature  is 
constant  or  not. 

Note  that  step  2  can  be  implemented  either  by  an  off-line  least-squares  fit  to  the  ICAO 
atmosphere  with  an  appropriate  model  or  with  a  table  and  an  appropriate  interpolation 
algorithm.  The  density  of  the  table  can  be  traded  off  against  accuracy  according  to  user 
needs.  A  power  series  expansion  good  to  *  0. 2  mb  when  using  36 -bit  floating  point 
is  as  follows: 

13 

»-  E 

i=  l 

A 

where  P  is  in  millibars, 


At  (0.4H)1"1  (4-7) 

H  is  in  kilofeet,  and  the  coefficients  Aj  are  given  in  Table  4-1. 
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An  accurate  interpolation  baais  is  as  follows: 

H~hk 

Vi  -  \ 

(4-8) 

where  and  h^+1  are  the  ICAO  table  entry  altitudes  just  smaller  and  just  larger  than  H 
respectively.  The  corresponding  pressures  in  the  table  are  pk  and  Plc+1* 

TABLE  4-1 

COEFFICIENTS  OF  ICAO  PRESSURE  MODEL 


1 

1.0133020  x  103 

2 

-9.2256754  x  101 

3 

4.6910444 

4 

-1.0669557 

5 

3.9830776  x  lO"1 

-2 

6 

-9.3698400  *  10 

-2 

7 

1.4027070  x  10 

6 

-1.3812160  x  10~3 

-5 

9 

9.0493196  x  10  J 

10 

-3.9030229  x  10'6 

-7 

11 

1.0639916  x  10 

12 

-1.6613529  x  10'9 

13 

1.1322430  x  10**11 

By  way  of  validation  as  well  as  illustration  of  the  method  consider  the  following 
example:  Figure  4-2  shows  actual  height  versus  range  detections  for  both  the  GE-692  radar 
and  the  IFF  Mode-C  link.  While  it  is  difficult  for  the  reader  to  unambiguously  follow  the 
track  for  each  sensor  on  the  two  inbound  ( I2  and  Ig)  and  two  outbound  legs  (0  g  and  0  4), 

the  following  table  provides  recorded  data  for  three  cases  of  Interest: 
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Figure  4-2.  Height  vs  Range  for  30  August  1979 


TABLE  4-2 


THREE  CASES  FOR  RADIOSONDE  CORRECTION 


CASE 

b 

b 

Ho  -  Hr 
%  c 

1 

27410 

26170 

+1240 

2 

47531 

46295 

+1236 

3 

7353 

6985 

+  368 

If  one  accepts  McUe-C  data  without  correction,  one  finds  quite  large  radar  bias  errors,  i.e. , 
H^-  Hc.  Table  4-3  summarizes  the  radiosonde  data  taken  that  day  at  two  surrounding  sites 
for  the  altitudes  of  interest. 


TABLE  4-3 

RADIOSONDE  DATA  FOR  AUGUST  30,  1979 


0  j 

ALBANY 

AVERAGE  1 

P 

H 

warn 

H 

■H 

H 

Tv 

(mb) 

(kit) 

1M 

(k ft) 

K3H 

(kft) 

(°K) 

300 

6.664 

282.44 

6.667 

285.17 

6.665 

283.81 

750 

i _ 

8.432 

279.45 

8.426 

281.58 

8.429 

280.52 

1 

2';.  645 

251.5 

24.568 

SOI 

ESSI 

252.25 

EES 

27.688 

243.5 

27.757 

ESS i 

BOB 

244.30 

|  150 

ESSI 

ESI 

42.814 

rmm 

46.039 

211.60 

j  125 

EH 

Di 

45.967 

49.732 

210. 6o] 

From  Equation  (4-7)  the  inferred  ICAO  pressure  at  thy  Kode-C  altitude,  H^,  is  given  in 
Table  4-4.  Usi&~  Equation  (4-5)  to  interpolate  from  Table  4-3,  we  get  ,  the  corrected 
Mode-C  height  as  given  in  Table  4-4.  The  radar  error  H^  -  H^,  has  been  reduoed  by 

TABLE  4-4 

THREE  CASES  AFTER  RADIOSONDE  CORRECTION 


CASE 

HC 
.  ffh 

pC 

(mb\ 

HC 

(tt) 

hr-h6 

tn\ 

1 

26170 

357.74 

27277 

+133 

2 

46295 

139.13 

47638 

-107 

3 

6985 

782.38 

7271 

+  82 

a  factor  of  10  at  the  two  higher  altitudes  (cases  1  and  2)  and  by  a  factor  of  4  for  the  third 
case.  While  not  conclusive,  this  does  provide  compelling  impetus  to  proceed. 


4.3  RESIDUAL  STANDARD  DEVIATION  OF  MODE-C  HEIGHT  BIAS  ERROR 
AFTER  RADIOSONDE  CORRECTION 

It  turns  out  that  the  radiosonde  data  table  can  be  accurately  interpolated  without 
reference  to  the  virtual  temperature  because  virtual  temperature  was  used  to  compute  the 
hydrostatic  altitudes  at  given  pressure  levels  in  the  first  place.  An  aocurate  interpolation 
formula  which  yields  comparable  results  for  the  three  cases  considered  in  the  previous 
section  ia, 


Hj  +  (H 


i  + 


-Hi> 


A 


in(Pnf/Pj> 

^  n  (pi+l/pi) 


(4-9) 


where 

H°.  is  the  hydrostatically  interpolated  Mode-C  height  (kft)  for  the 

m  data  point  and  i  flight  leg, 

P.  is  tho  inferred  altimeter  pressure  (mb),  logarithmically 

interpolated  from  the  ICAO  standard  atmosphere  at  altitude  h^  , 
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Pj  and  p. 


i+1 


H,  and  Hi+1 


are  the  bracketing  pressures  (mb)  in  the  radiosonde  data, 

are  the  hydrostatic  altitudes  (kft)  corresponding  to  pt  and  p^+1,  and 

is  the  received  Mode-C  height  (kft). 


The  inferred  altimeter  pressure  Pni  is  extracted  from  the  ICAO  standard  atmosphere  at 
altitude  h^  from, 


/s 


(4-10) 


where 


h^  and  h^+1  are  the  bracketing  altitudes  (kft)  in  the  ICAO  table,  and 

P,  and  P,  -  are  the  ICAO  standard  pressures  (mb)  which  correspond 

tol^andh^. 


The  objective  of  this  section  Is  to  characterize  the  residual  error  in  the  Mode-C 
height  after  corrections  derived  from  applicable  radiosonde  data,  as  prescribed  here  and 
in  pa.".  4. 2.  The  principal  sources  of  error  in  the  correction  procedure  derive  from: 

1.  The  error  derived  from  the  measurement  and  use  of  radiosonde  data,  and 

2.  Tiio  eu’or  in  retrieving  the  true  pressure  as  perceived  by  the  altimeter 
barometer. 

From  Equation  (4—7)  we  can  relate  errors  in  radiometer  pressures,  pj  and  p^+^, 

and  the  inferred  ICAO  standard  pressure,  P  .  ,  to  the  corrected  Mode-C  height,  H° . 

nx  m  • 

Thus, 


s»n°<  1 

r  Hi+ii 

Hnf 

<SPt 

rvw 

A 

3P  . 
nl 

r.  5.] 

ft±i 

»ni  in/Pi±l\ 

A  o 

<1 

m  mi 

M 

A 

pnl 

*  ^  c 

pi+l 

(4-11) 
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Now  we  can  relate  the  error  in  the  inferred  ICAO  standard  pressure  to  measurement  and 
encoder  errors  from  Equation  (4-10). 


5henc  f Pk+1 


-  in  I - 

\+l~\  \Pl 


(4-12) 


where 


<5h  is  the  encoder  error  (  o.  *  65  ft),  and 

sue  n _ 

enc 

6  Pm  is  the  pressure  measurement  error 

(ap  =  3  mb), 
m 

Lenhard  in  reference  3  investigated  the  uncertainties  associated  with  the  radiosonde 
procedure.  Unfortunately  he  only  characterized  the  standard  deviation  component  of  error 
and  what  follows  Is  similarly  limited.  Here  we  regard  his  results  as  impacting  the  standard 
deviation  of  the  residual  error  in  the  Mode-C  height  after  correction  by  radiosonde  means. 
His  results  are  summarized  in  Figures  4-3  and  4-4  which  plot  the  standard  deviation  of 
radiosonde  computed  height  versus  pressure  and  the  standard  deviation  of  radiosonde 
computed  pressure  versus  height  respectively.  From  Figure  4-4  we  get  dp  directly  versus 
altitude.  Presuming  independence  between  all  sources  of  uncertainty,  we  combine  Equations 
(4-11)  and  (4-12)  and  take  the  expectation  value.  Thus, 


<i) _ i _  [A . 

t  ( Hi  -  Hi+i\2  [Yi“ 

\  Ki  )  LVv 


ti  \2  .  w  \*  /  „  \2-  /  \  2 

K'J  v  W ) 

.EJs+L  \  /a  \ 2  *1 

-*-)  4  ,+wA  .  «-> 

,-h.  /  enc  \Pni  /  J 
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Figure  4-3.  Standard  Deviation  of  Radiometer  Computed  Height 
Error  vs  Pressure  -  after  Lenhard  (Table  4-1) 


Figure  4-4.  Standard  Deviation  of  Radiosonde  Computed  Pressure  at  a 
Specified  Altitude  (Inverse  Hydrostatic  Problem)  -  after 
Lenhard  (Table  4-2) 


This  error  is  plotted  In  Figure  4-5  as  the  upper  curve.  The  second  curve  down  shows  the 
effect  of  encoding  error  while  the  third  shows  the  effect  of  the  radiosonde  induced  errors 
alone.  The  primary  component  of  tie  error  plotted  in  Figure  4-5  is  due  to  the  extraction 
of  ICAO  standard  pressure  from  the  Mode*C  height  and  the  attendant  altimeter  pressure 
measurement  error.  As  such,  the  greatest  advantage  of  the  radiosonde  correction  is  at 
altitudes  below  about  45  kft  as  can  be  seen  by  comparing  the  plots  of  Figures  3-5  and  4~G. 
Fortunately,  this  is  the  important  region  for  flight  tests.  Nevertheless  the  residual  error 
after  correction  at  45  kft  still  has  a  standard  deviation  of  about  450  ft. 


Figure  4-5.  Standard  Deviation  of  Combined  Altimeter  and 
Radiosonde  Height  Error  vs  Height 
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SECTION  V 


ALTERNATIVE  APPROACHES  TO  HEIGHT  BIAS  TESTING 

In  the  previous  section,  we  formed  the  basis  of  an  accurate  height  reference  by 
introducing  a  method  of  correction  to  the  Mode-C  height.  Here  we  expose  two  methods  for 
testing  radar  height  bias  using  this  corrected  Mode-C  height  reference. 

5.1  FIXED  HEIGHT  BIAS  TESTING 

A  hypothesis  is  advanced  that  all  radar  bias  errors  are  random  in  the  long  run  and 
zero  mean.  The  test  is  designed  to  sample  this  composite  long-term  distribution  by 
averaging  the  difference  between  the  radar  (primary)  and  the  IFF  Mode-C  (secondary)  over 
many  flights  of  common  purpose,  e.  g.  in -the -clear  at  40  kft.  If  radiosonde  data  is  avail¬ 
able,  it  is  used  to  compensate  and  correct  the  secondary  reference  altitude  on  each  flight 
leg.  Pass/fail  thresholds  are  computed  from  the  following  constraints: 

1.  The  long-term  average  should  be  zero  in  the  absence  of  fixed  bias  errors, 

2.  All  bias  tests  flights  are  far  apart  enough  in  time  so  that  independent  bias 
samples  are  taken  from  radar  equipment  and  atmospheric  sources. 

3.  Producer's  and  buyer's  risks  are  bounded  by  accommodating  the 
appropriate  sampling  uncertainty. 

4.  If  no  radiosonde  data  is  available,  the  expected  variance  of  the  altimeter 
error  is  Included  in  the  calculations  of  pass/fall  thresholds. 

5.2  TOTAL  MEAN-SQUARE -HEIGHT  ERROR  TESTING 

The  first  approach  implies  separate  pass/fail  criteria  and  thresholds  for  height 
bias  from  those  used  to  test  the  variance  of  height  accuracy.  The  second  test  approach 
attempts  to  combine  the  two  error  sources  (i.  e.  bias  and  short-term  random)  into  a  test 
of  the  mean-square  absolute  height  error. 

The  squared  difference  of  primary  and  secondary  heights  Is  averaged  on  special 
flights  designated  for  bias  tests.  Altimeter  corrections  are  made  on  eaoh  flight  leg  in 
range  intervals  where  radiosonde  data  is  available.  The  sample  varianoe  is  adjusted  for 
the  expected  variance  of  atmospheric  refraction  biases  (if  refraction  bias  is  not  a  respon¬ 
sibility  of  the  radar  producer). 


The  same  constraints  hold  for  each  test  type  and  pass/fail  thresholds  are  oomputed 
from  knowledge  of  expected  theoretical  performance  with  allowances  for  the  limited  sample 
size.  In  order  to  estimate  the  tolerance  required  to  accommodate  the  uncertainty  associated 
with  a  finite  sample  size,  it  is  necessary  to  understand  the  statistics  of  the  sample  estimate 
of  both  the  fixed  height  bias  error  and  the  mean-square  total  height  error. 

5.3  ERROR  ANALYSIS  FOR  THE  FIXED  HEIGHT  BUS  TEST 

In  this  section  we  will  define  the  error  sources  and  introduce  correlation  between 
samples.  Inpa; .  5. 4  we  will  then  derive  expressions  for  producer  and  buyer  rlskB  for  test 
planning.  In  par.  5. 5  we  will  consider  this  test  in  the  context  of  an  overall  flight  teBt  where 
ether  performance  measures  are  concurrently  tested  so  as  to  introduce  the  concept  of  risk 
allocation.  Numerical  examples  to  tie  all  these  ideas  together  are  included  in  par.  5. 6. 

5.3.1  DEFINITIONS 

As  a  prelude  to  the  analysis,  consider  the  definitions  for  a  given  flight  test  type 
(e.  g.  in -the -clear  at  40  kft)  and  a  given  flight  path  interval.  For  the  I  ^  flight  leg  and  n^ 
data  point,  the  radar  (primary)  errors  are, 

a .  the  total  short-term  random  error  component  of  elevation  angle 

(jitter  and  thermal).  It  is  independent  from  point-to-point  within 
the  estimation  interval  and  is  generally  considered  as  zero-mean 
Gaussian,  whose  variance  c r2  may  vary  systematically  with  range 
and  elevation  angle,  ° 

e,  the  total  equipment  (electronic,  mechanical,  and  software) 

elevation  bias  error  which  is  random  and  independent  from 
leg-to-leg.  Its  parent  distribution  is  assumed  to  be  zero- 
mean  Gaussian  with  variance  ffc2  ,  which  may  vary 
systematically  with  range  and  elevation  angle, 

A.  the  total  atmospheric  elevation  bias  error  which  is  random  and 

independent  from  leg-to-lej'.  Its  parent  distribution  (after 
refraction  corrections)  is  assumed  to  be  zero-mean  Gaussian 
with  variance  that  varies  systematically  with  elevation  angle, 

the  total  zero-mean  Gaussian  bias  error  which  is  hi  dependent 
from  leg-to-leg.  It  has  variance  erg  ■  , 

c  constant  bias  error  with  zero  variance; 

<5  HyAR  Mode-C  error  without  radiosonde  correction  in  the  absenoe 
of  altimeter  error,  assumed  independent  from  leg-to-leg 
and  zero-mean  Gaussian  with  variance  , 

c  the  height  bias  estimate, 


H^q  radar  height  measurement  for  the  point  and  i  ^  leg, 

til 

H?  uncorrected.  Mode -C  height  measurement  for  the  n 

1 11  point  and  iTh  leg, 

H  mean-range  in  j/ath  intervals, 

5Ha  r  rp  Mode-C  error  without  radiosonde  correction  in  the  absence  of 
AL|1  atmospheric  pressure  variability,  assumed  Independent  from 

point-to-point  and  zero-mean  Gaussian  with  variance  t 

ALT 

fiHRAD  Mode-C  height  error  after  radiosonde  correction  for  zero 
KAU  altimeter  error  which  are  assumed  Independent  from  leg-to-leg 
and  Is  zero-moan  Gaussian  with  variance  ,  and 

<5Hjw  Mode-C  height  error  after  radiosonde  correction  in  the  absence 
of  radiosonde  error  which  is  independent  from  point-to-point 
and  zero-mean  Gaussian  with  variance  o|  ^  . 

5. 3. 2  VARIANCE  OF  THE  ESTIMATE 

The  height  bias  estimate  after  averaging  over  the  N  data  points  in  the  path 
interval  and  L  independent  flight  legs  Is 


L  N  v 

E  E  {Eln  -  Hi°n) 


l  =1  n=l 


where 


H?n 


(“in  +  3i)  +c+  (4HVA»)1  +  (5HALt), 


if  the  altimeter  is  not  corrected  by  radiosonde,  or 


*fn  "  Hfn  =  *  (°fn  +  Pij  +  c  +  (^Ra), Q  4 


if  it  is.  Thus, 


c  «  c  + 


T5T  ^  E  H“in+  (aHRA)in  +  (shalt)J; 

i  E  [ &  si  +  (5Hvar).  +  (4HRAd). 
i=i  L 


(5-1) 


(5-4) 
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where  <5H^LT  and  are  Ignored  If  radiosonde  is  used  and  the  errors  and 

<5Hj^ad  are  ignored  otherwise.  Note  that  the  estimate  of  o  is  unbiased  in  that 


E  (6)=  c 
because 

E  («)-  E  (/3)  =  E(6HVAR)  =  0 
E(jHALT>sE(iHRAD>  “E<6HRA>“° 
by  definition. 

The  variance  of  the  estimate  is, 


VARc  =  |r2VAR«  +  VARfdHj^J  +  VAR  (6Halt)  I 

L  1  (6“ 
+  |"r2  VAR/3  +  VAR (^AD)  +  VAR  (aHvar)J 

where  all  error  sources  are  assumed  to  be  uncorrelated  by  type.  The  standard  deviation 
of  the  unbiased  estimate  c  is 


JL  r2q.2 

LN  LR  a 


2  +  +  or? 


°RA  +  °ALT 


“f  n  ■ 

.  +  T  [**  4  *  "Iad  +  °vaeJ  .  <5-7> 


The  estimate  c  is  both  unbiased  and  consistent  because 


E(c-c)  *  0 


(6-8) 


lim  <r~  *  0. 
c 

L-*  00 


(5-9) 


5.3.3  SAMPLE-TO-SAMPLE  CORRELATION 


Equation  (5-7)  is  correct  only  if  each  leg  is  flown  far  enough  apart  in  time  to  inde¬ 
pendently  sample  all  the  bias  errors.  To  gain  some  insight  into  the  penalty  to  be  paid  for 
failure  to  do  this  consider  the  following  informal  analysis: 

Let: 

T  *  flight  time  of  one  leg, 

AT  time  between  legs, 

r  total  bias  correlation  time, 

L  total  number  of  legs, 

Leff  effective  number  of  independent  legs,  and 
p  correlation  coefficient  of  bias  errors. 

A  reasonable  model  for  versus  bias  correlation  can  be  constructed  from  the 
following  boundary  conditions: 

1.  For  zero  correlation,  -  L.  d  Lef(/dp=0. 

2.  For  unity  correlation,  Leff  -  »• 

Without  taking  this  quest  rigorously,  we  can  satisfy  these  criteria  with  a  simple  cosine 
function,  i.e. 

Loff-  -§=*-  Cos  1TP  +  jf-  (5-10) 

which  is  plotted  in  Figure  5-1  as  the  ratio  of  L^/L  with  L  as  a  parameter.  Th’s  plot 
implies  that  a  correlation  as  low  as  0.  25  is  requred  to  achieve  90%  independence  between 
legs. 


A  good  model  for  correlation  as  a  function  of  the  total  flight  cycle  time  T  +  AT  is 
the  Markov  process.  If  the  current  bias  error  only  depends  upon  the  previous  error  (1.  e. 
the  value  on  the  last  flight  leg)  then  we  are  describing  a  1st  order  Markov  process  with 


P 


T+AT 

*  J  * 


(6-11) 


Figure  5-1.  Effective  Number  of  Independent  Legs  vs  Correlation 
With  Number  of  Legs  as  a  Parameter 


If  it  depends  upon  the  last  two  errors  (legs)  then  we  are  describing  a  2nd  order  Markov 
process  with 


(5-12) 


These  are  plotted  in  Figure  5-2  and  combined  with  Equation  (5-10)  in  Figure  5-3.  Here  we 
see  90%  independence  can  be  achieved  with  a  total  flight  cycle  time  of  about  1,5  t,  I  have 
not  uncovered  definitive  information  on  t  but  experience  with  the  GE-592  implies  2-8  hours. 
So  to  insure  independence,  legs  for  the  bias  test  should  be  flown  no  more  than  twice  a  day 
(i.  e.  about  12  hours  apart). 

To  be  precise  we  should  modify  the  expression  for  estimation  error,  i.  e.f  Equation 
(5-7),  to  include  leg-dependence.  The  first  term  includes  errors  which  are  independent  from 
point-to-point  so  that  =  L.  But  the  last  two  terms  are  bias  errors  so  that 
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5.  3.  4  ERROR  MODELS 


To  gauge  the  value  of  these  ideas,  we  need  representative  models  for  and 
for  a  typical  air  defense  radar  (not  necessarily  the  GE-592)  without  special  features  (e.  g. 
multipath-defeating  angle  measurement  techniques,  SNR  weighting  for  multiple-hits,  etc. ). 
For  random  errors, 
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=  a 


CL 


+  K 


with 


0.65  mrad;  R  £  R 


0. 9  mrad;  R  >  R 


(5-14) 
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2ND  ORDER  MARKOV  PROCESS 

p  =  e"2*146  <T+AT)/r  +  2.14s  (  T  +  aT)/t) 


■■mr1 


1ST  ORDER  MARKOV  PROCESS 
P  -  exp  (-(T+aT)/t  ) 
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Figure  5-2.  Correlation  as  a  Function  of  the  Ratio  of  Flight 
Cycle  Time  to  Bias  Correlation  Period 


The  error  variance  °VAR  is  identified  as  the  Mode-C  height  error  variance  without 


radiosonde  correction  and  in  the  absence  of  altimeter  error.  The  sources  of  uncertainty 
are  the  temporal  and  spatial  variability  of  pressure  at  the  surface  and  at  altitude.  It  seems 
reasonable  to  assume  that  this  error  source  is  independent  only  from  leg-to-leg.  Again 
from  Equation  (3-6), 


where, 


0. 032  pQ  , 


(5-22) 

(5-23) 


and 


a  =  204.8  kft, 
b  =  (103. 33  kft)2, 
c  =  36.  922  mb/kft, 
d  =  (269. 503) 2  mb2 /kft,  and 

P0«  1013. 25  mb. 

o 

The  error  variance  is  the  variance  of  the  Mode-C  height  error  after  radio¬ 

sonde  correction  in  the  absence  of  altimeter  error.  Thus  from  Equation  (4-1), 


where  cr^  Is  the  variance  of  the  smoothed  radiosonde  pressure  at  a  given  hydrostatic 
altitude  as  given  in  Figure  4-4.  Note  that  this  error  is  independent  only  from  leg-to-leg. 

The  error  variance  cr2  is  the  variance  of  the  Mode-C  height  error  after  radiosonde 

P 

correction  in  the  absence  of  radiosonde  error.  This  error  is  independent  from  data  point 
to  data  point  and  also  comes  from  Equation  (4-1)  as, 


HrHi+i 


i  n  < 


i+1 


(5-25) 


The  parameters  which  materially  affect  the  uncertainty  of  the  fixed  height  bias 
estimate  are: 

1.  The  ratio  of  the  leg  cycle  time  to  the  leg-to-leg  correlation  time, 
i.e. ,  (T+AT)/t. 

2.  The  number  of  flight  legs,  L  , 

3.  The  flight  altitude,  H. 

In  order  to  expose  these  sensitivities,  constant  altitude  radial  flight  legs  were  simulated 
with  a  target  speed  of  about  500  knots  and  a  data  interval  of  12  seconds.  This  gives  about 
7  data  points  per  bias  estimate,  per  leg.  In  Figure  5-4,  ten  legs  were  assumed  with  the 
aircraft  at  40  kft.  Three  ratios  of  (T+  AT)/r  were  selected,  1.  e.  0. 1,  0. 6,  and  1. 6 
and  the  second-order  Markov  model  was  used  for  correlation.  As  it  happens  a  ratio  of 
1. 5  is  almost  equivalent  to  an  infinite  number  of  oor r e latloo -time s  between  legs,  or  a 
leg-to-leg  correlation  of  xero.  Note  the  dramatic  effect  of  radiosonde  correction  for  all 
ranges  below  about  120  nmi.  These  curves  say  that  if  one  is  only  concerned  about  fixed 
height  biases  of  a  few  thousand  feet,  then  it  probably  isn't  neoessary  to  use  radiosonde 
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Figure  5-4.  Standard  Deviation  of  Fixed  Height  Blaa  Estimate 
with  and  without  Radiosonde  Correction 
-  Effect  of  Correlation  from  Leg  to  Leg 


"■  ri»W>CT  VP1 TI 


so  long  as  T  +  AT  a  1,5  t.  On  the  otherhand,  If  fixed  height  biases  of  a  few  hundred 
feet  can  be  worrisome  then  radiosonde  correction  Is  essential. 

The  effect  of  the  number  of  legs  Is  shown  In  Figure  5-5  with  T  +  A  T  a  1. 5  r, 

N=7,  and  H*40  kft.  Here  we  can  see  the  value  in  using  10  legs  but  probably  not  for  30, 
particularly  with  radiosonde  correction.  (This  question  may  require  revisiting  when  we 
consider  the  statistical  risk). 

Finally  consider  the  effect  of  aircraft  altitude.  Results  for  altitudes  of  10,  20,  30, 
and  40  kft  are  plotted  in  Figure  5-8  where  N=7,  u*10,  and  (T  +  AT)  a  1. 5  r.  The  strange 
resonance-like  behavior  at  10  and  20  kft  derive  from  the  fact  that  at  these  peaks  the  eleva¬ 
tion  angle  is  zero  and  the  uncertainty  from  atmospheric  refraction  is  dominant.  This 
effect  is  so  dominant  that  radiosonde  correction  is  totally  useless.  Of  course  for  a  sea 
level  site  the  far  side  of  these  peaks  cannot  be  seen  by  the  radar  as  they  correspond  to 
negative  elevation  angles.  Nevertheless  these  curves  do  provide  at  least  qualitative  effects 
for  raised  sites  where  negative  angles  are  possible^ 

5.4  STATISTICAL  RISKS 

Since  the  error  dletributlen  of  3  -  o  I,  hero-mean  Oeuealan  with  variance  <r§  (Eqint- 
tion  (5-13)),  pass/fail  thresholds  are  easily  formulated  in  terms  of  risks.  The  producer's 
maximum  risk  is  the  probability  of  failing  the  test,  i.  e.  |  c  |  >  |cp/p  |  where  |Cp  |  is 

the  pass/fail  threshold,  given  that  the  true  fixed  bias  error  (in  tho  absenoe  of  sampling  un¬ 
certainty)  is  zero,  i.  e.  false  failure.  Thus, 

RS  “  Pr  [l S  I  >  t  Cp/F  1/  !  0  I  "  °]  •  (5-26) 

The  buyer's  maximum  risk  is  the  probability  of  passing  the  test,  1.  e.  j  8  |a»  |  Cpyp  |  , 

given  that  the  true  fixed  bias  error  ]  c  |  is  his  maximum  tolerated  fixed  bias  error  )CB  | , 
i.  e. ,  false  passing.  Thus, 

RS  -  Pr  [l®  Is  I  CP/F  |/  1°  I  -  |CB  I  ]  (S-27) 

These  risks  are  illustrated  in  Figure  5-7.  Evaluated  in  terms  of  error  functions. 
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Figure  5-5.  Standard  Deviation  of  Fixed  Height  Bias  Estimate 
with  and  without  Radiosonde  Correction 
-  Effect  of  Number  of  Legs 
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Figure  5-7,  Schematic  Distributions  of  Estimation  Error  for  Producer 
and  Buyer  Showing  Risks  as  Shaded  Areas 
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These  risks  are  plotted  in  Figures  5-8  and  5-9,  respectively,  for  typical  ranges  of  the 
parameters  |Cp^F  j  /a^  and  }  C0  |  /  <Tg.  Note  that  once  a  producer's  risk  Is  allocated 
from  an  overall  risk  level  down  to  the  marginal  risk  for  the  height  bias  test,  *hi»t  the  pass/ 
fall  threshold  onljr,  depends  on  the  uncertainty  of  the  estimate,  1„  e.  as  measured  by  <rg,  the 
quantity  evaluated  In  the  previous  section.  The  buyer's  position  Is  not  symmetrical  with  the 
producer's  position  however  In  that  his  least  acceptable  performance  (as  represented  by  the 

buyer's  Umlt  JCfi  J )  depends  on  the  pass/fall  threshold  as  well  as  o>  ,  for  a  given 
buyer's  risk.  0 
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Figure  5-9.  Buyer's  Risk  vs  Normalized  Maximum  Tolerable 
Fixed  Biaa  Error  with  Normalized  Pass/Fall 
Threshold  as  a  Parameter 
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Before  fixing  ideas  further  with  numerioal  samples,  it  is  helpful  to  understand  the 
risk  allocation  process  so  that  we  can  assess  realistic  risk  levels. 

5.5  OVERALL  TEST  STRUCTURE 

Before  we  can  go  back  to  the  application  where  we  quantified  typical  values  for  the 
uncertainty  in  the  bias  estimate  and  to  apply  these  values  to  help  gauge  pass/fail  thresholds, 
we  must  first  assess  the  levels  of  reasonable  risk.  A  mutually  acceptable  range  of 
overall  acceptance  test  risks  appears  to  be  between  5  and  10%.  That  is,  the  risk  of  the 
producer  falsely  failing  or  falsely  passing  the  controlled-flight  acceptance  test  merely 
because  of  inadequate  sample  sizes  is  between  5  and  10%.  The  bias  test  is  only  one  compo¬ 
nent  of  the  overall  test.  In  fact,  as  many  as  14  Independent  performance  measures  can  be 
tested,  i.  e. , 

1.  Probability  of  detection 

2.  Absolute  range  accuracy 

3.  Absolute  height  accuracy 

4.  Absolute  azimuth  accuracy 

5.  95th  percentile  absolute  height  accuracy 

6, 7.  Relative  height  accuracy  between  closely  spaced  and  distantly 
spaced  targets 

8, 9.  95th  percentile  relative  height  accuracy  between  closely  spaced 

and  distantly  spaced  targets 

10, 11.  Range  resolution  between  targets  of  equal  and  unequal  radar  cross 
section  (RCS) 

12, 13.  Azimuth  resolution  between  targets  of  equal  or  unequal  radar 
cross  section  (RCS) 

14.  Height  (or  elevation)  bias  test 

where  absolute  accuracy  refers  to  a  single  aircraft  and  can  include  bias  in  an  rms  sense  or 
just  involve  the  variance.  Relative  accuracy  implies  two  aircraft.  There  can  be  as  many 
as  four  separate  altitudes  or  types  of  environmental  conditions  tested  and  separate  pass/fail 
thresholds  established  for  as  many  as  20-range  Intervals  per  radial  leg.  Thus,  there  can 
be  hundreds  of  separate  pass/fail  thresholds  for  an  acceptance  test  with  an  overall  risk 
between  5  and  10%. 


5. 5. 1  PASS/FAIL  CRITERIA 

What  must  tbe  allocated  risk  be  on  the  level  of  a  single  range  interval  (where  pass/fail 
thresholds  are  applied)  ?  This  depends  upon  the  criterion  for  passing  the  overall  test.  If 
the  producer  must  pass  each  and  every  subtest  to  pass  the  acceptance  test,  then  the  pro¬ 
ducer's  risk  will  add  up  rapidly  and  his  risk  per  subtest  must  be  very  small  indeed.  If  the 
allocated  risk  per  lowest  level  suptest  rp  is  equal  for  all  subtests,  then  for  an  overall  risk 
level  of  Rp, 


Rp  «  1  -(1  -  rp)  Ns 


(5-31) 


where  Ng  is  the  number  of  subtests.  Similarly  for  the  buyer,  there  are  so  many  ways  for 
the  producer  to  fail  that  the  buyer's  risk  per  lowest-level  Buhtest  will  reduce  rapidly.  If 
the  overall  risks  to  buyer  and  producer  are  to  be  about  equal  (seems  reasonable)  then  the 
buyer's  subtest  risk  rB  must  be  fairly  large.  In  fact 


n  Ns 
=  (rB) 


(5-32) 


These  equations  are  easily  solved  for  and  rB  and  are  plotted  in  Figures  5-10  5-11  as 

the  curve  marked  M=0,  where  M  is  the  to*al  number  of  subtest  failures  permitted  for  the 
producer  to  pass  the  acceptance  test.  Note  that  for  10  subtests  the  producer's  risk  per  sub¬ 
test  must  be  as  low  as  1%  to  achlave  an  overall  level  of  10%.  Moreover,  for  Ns  «=  100, 
rp  »  0. 1%.  Similarly,  for  10  subtests  the  buyer  must  endure  a  riu.1'  per  subtest  of  almost 
80%.  And  when  Nfl  =  100,  rB  ~  98%.  Buyers  often  chafe  at  this  notion  (although  perfectly 
true)  and  seek  to  relieve  this  pressure  while  simultaneously  raising  the  producer's  sub¬ 
test  risk  by  permitting  the  produoer  to  fail  a  few  subtests  while  still  passing  the  acceptance 
test  itself.  Quantitatively, 
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The  remaining  curves  in  Figures  5-10  and  5-11  show  the  effects  of  permitting  1  to  10  total 
subtest  failures.  An  interesting  way  to  gauge  the  effectiveness  of  this  prooedure  is  to  con¬ 
sider  an  exercise  whereby  the  producer's  risk  per  subtest  remains  at  the  overall  level 
(say  10%).  The  first  column  of  the  following  table  shows  the  required  number  of  tests  for 
this  to  occur  for  M  between  0  and  10.  The  third  and  fourth  columns  give  the  corresponding 
producer's  and  buyer's  subtest  risks.  Compare  these  to  the  last  two  columns  which  must  be 
obtained  for  no  permitted  failures.  The  contrast  is  marked.  As  we'll  see  shortly,  the 
corresponding  pass/fail  thresholds  and  the  buyer's  maximum-tolerated  fixed  bias  level  will 
also  be  affected. 
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10 

77 

0 

0. 35 

93.8 

70 

10 

10 

80 

0 

0. 15 

96.5 

5.5.2  PASS/FAIL  THRESHOLDS 

Now  we  can  combine  the  estimation  uncertainty  c.  and  pass/fail  criteria  into  pass/ 
fail  thresholds  and  buyer's  limits.  Figure  5-12  is  a  plot  vZ  ihe  pass/fail  threshold  versus 
range.  Here  we  have  selected  the  case  'where  the  flight  cycle  time  Is  long  enough  with 
respect  to  the  bias  error  correlation  period,  1.  e. ,  T  +  AT  a  1. 5  r,  that  the  legs  are  effec¬ 
tively  independent.  Moreover  we  consider  the  case  for  10  legs  with  and  without  correction 
to  the  Mode-C  height  and  three  levels  of  produce's  subtest  risk,  i.  e. ,  0. 1%,  1%,  and  10%. 
From  Figure  5-12  we  can  see  a  dramatic  difference  St  close  ranges  with  and  without  radio¬ 
sonde  correction.  For  example,  at  60nmi  aadr^«0. 1%,  f  Cp/jp  |  *  1400  ft  without 
correction  and  only  340  ft  with  compensation,  almost  a  factor  of  4.  At  200  nmi,  however, 
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Figure  5-12.  Pass/Fail  Threshold  vs  Range  for  Fixed  Height  Bias 
Error  at  Producer's  Risks  of  0. 1,  1,  and  10% 
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the  estimation  error  is  dominated  by  thermal  errors  which  cannot  be  adequately  reduced 
even  with  7  data  points  per  range  internal  and  10  legs  (intervals  per  subtest)  for  a  total  of 
70  samples.  Thus,  for  r  -  0. 1%  |Cp/F  f  ■  2850  ft  without  compensation  and  still  as  high 
as  2500  ft  even  after  altimeter  correction.  The  value  of  radiosonde  corrections  is  then 
questionable  for  ranges  beyond  about  120  nmi. 

The  question  of  the  maximum  fixed  bias  the  buyer  must  tolerate  to  accommodate 
sampling  uncertainty  is  addressed  next.  Consider  two  contrasting  criteria  for  passing: 

1.  M=0  (no  permitted  failures) 

P  B 

For  R  =R  *  10%,  from  Figures  5-10  and  5-11 

rP  *  0. 15%;  rB  -  96. 5%. 

From  Figures  5-8  and  5-9, 

ICP/F  |  ■  3.20  ffg  ;  \  CB  !o  *  1.4  Og  , 

2.  M«10  (10  out  of  70  permitted  failures) 

For  RP=RB  «  10%,  from  Figures  5-10  and  6-11 

r^j  e  10%;  rBQ  »  80%. 

From  Figures  5-8  and  5-9, 

! cp/f I  “  l-6  a8  ‘  I°b!io  "  °-7V 

Here  we  see  an  interesting  phenomenon.  The  strategy  of  permitting  some  subtest  failures 
does  four  good  things  for  the  buyev.  First,  it  forces  the  producer's  subtest  risk  up  from 
0. 15%  to  10%.  Second,  this  forces  the  pass/fall  threshold  down  from  3. 2  o~  to  1. 6  erg  , 
a  factor  of  2.  Third,  the  buyer's  subtest  risk  reduces  from  96. 5%  down  to  80%.  Finally, 
the  buyer's  limit  reduces  from  1, 4  a*  for  no  permitted  failures  to  0. 7  Og  for  10  out  of  70 
permitted  failures,  also  a  factor  of  2.  The  buyer  will  usually  have  free  choioe  here,  tail¬ 
oring  his  strategy  to  his  particular  interests  and  concerns.  From  the  producer's  standpoint, 
his  total  risk  is  being  guaranteed  in  any  case  and  he  should  be  satisfied  to  give  the  buyer  this 
degree  of  freedom.  Figures  5-13  and  5-14  are  plots  of  |  Cp^p  j  and  |  CB  | ,  respectively, 
for  the  two  alternative  strategies  highlighted  here.  Once  again,  a  striking  conclusion  is  the 
limited  benefit  of  compensating  the  Mode-C  height  beyond  about  120  and  the  tremendous 
desirability  for  doing  it  at  oloser  ranges. 
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Figure  5-13.  Pass/FaJl  Thresholds  vs  Range  for  70  Subtests  with  Zero 

and  10  Permitted  Subtest  Failures  and  an  Overall  Risk  of  10% 


Figure  5-14.  Maximum  Tolerable  Buyer's  Limit  of  Fixed  Height 
vs  Range  for  70  Subtests  with  Zero  and  10  Permitted 
Failures  and  an  Overall  Risk  of  10% 
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5.5.3  TEST  CONDUCT 


Sometimes  only  one  radiosonde  will  be  available  in  the  vicinity  of  a  useful  radial 
flight  path  and  at  appropriate  times.  Experience  has  shown  that  reasonable  corrections  to 
the  altimeter  can  be  made  for  ±15  nml  about  the  radar  range  to  the  balloon.  Figure  5-15  is 
a  Range-Height-Angle  chart  which  assumes  that  the  radiosonde  balloon  ascended  in  the  70  to 
80  nmi  range  window.  Loosely  interpreting  this  to  a  window  of  applicable  altimeter  correc¬ 
tions  of  from  60  to  90  nmi,  we  can  see  that  four  flight  altitudes  can  cover  all  beams  from  0. 5° 
to  7°  in  elevation.  Moreover,  if  a  waveform  change  occurs  at  80  nmi,  these  tests  (at 
altitudes  of  10,  18,  30,  and  47  kft)  will  also  test  the  height  bias  in  both  short-range  and  long- 
range  beams. 

Because  we  require  no  more  than  two  bias  flights  a  day  (10-12  hours  apart)  at  a  given 
altitude,  a  sensible  approach  is  to  intersperse  the  bias  flights  with  all  other  flights.  For 
instance,  at  the  end  of  each  day  before  returning  to  base,  the  pilot  can  fly  a  single  leg  at 
each  of  the  four  altitudes  (spiraling  up  or  down  between  legs)  in  about  the  same  tima  it  takes 
to  fly  one  leg  from  0  to  200  nml.  Separate  statistics  are  accumulated  from  day-to-day  and 
separate  thresholds  applied  at  each  of  the  four  altitudes  (in  each  10  nmi  range  interval  or  in 
one  30  mni  interval)  at  the  end  of  the  allotted  number  of  legs  for  the  bias  test. 

6.6  ERROR  ANALYSIS  FOR  THE  MEAN-SQUARE  HEIGHT  ERROR  TEST 

In  par.  5. 3  through  5. 5,  we  analyzed  a  test  for  uncovering  a  fixed  bias  error,  the 
kind  which  results  from  an  error  in  calibration,  algorithmic  compensation,  or  software. 

The  buyer  is  primarily  interested  In  this  kind  of  error  and  already  accepts  the  budgeted 
equipment  bias  error,  which  is  random  from  day-to-day,  i.  e. ,  N(0,  a  ).  He  is  protected 
from  larger  than  promised  values  of  because  the  pass/fail  thresholds  are  calculated 
from  that  value.  The  second  test  type  accomplishes  similar  goals,  but  involves  a 
threshold  on  the  total  mean  square  height  error  in  the  30  »mt  window  of  applicable  radio¬ 
sonde  corrections.  In  both  oases,  the  variance  of  absolute  height  accuracy  is  estimated  and 
thresholded  separately  in  the  remaining  range  intervals.  In  section  we  develop  the 
statistics  of  the  total  mean  square  error  for  this  second  test  type  and  oompute  risks  and 
pass/fail  thresholds  to  exemplify  and  clarify  the  approach.  Mitch  of  what  we  have  already 
covered  remains  applicable. 
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Figure  5-15.  Typical  Elevation  Coverage  for  a  Height  Bias  Test 

with  a  Radiosonde  Ascent  Between  70  and  80  nmi  Range 


5.6.1  DEFINITIONS 


MSE  '  total  mean-square  error  estimate, 

w  fk  fk 

Hj  the  height  residual  for  the  n  data  point  and  i  leg, 

i.  e. ,  H*  -  H°,  where  the  altimeter  has  been 
corrected  for  temperature, 

ag  standard  deviation  of  the  height  residual, 

A,  a,  b,  v*  parameters  of  the  noncentral  chi-square  distribution, 

noncentral  chi-square  random  variable, 

standard  chi-square  random  variable  with 
v  degrees-of-freedom, 

f  (x)  probability  density  function  of  the  random  variable  x 

x  at  the  value  x,  and 

P  cumulative  probability  distribution  of  the  random 

x  variable  x. 


5.6.2  VARIANCE  OF  THE  ESTIMATE 

The  total  mean-square  error  is  the  average  squared  difference  between  the  primary 
and  secondary  height  predictions  with  an  adjustment  for  the  residual  root-mean-square  (rms) 
refraction  error,  <r^,  if  not  included  in  the  budget,  and  an  altimeter  presumed  to  be  error- 
free  (corrected  for  surface  pressure  and  by  radiosonde)41.  Thus 


A  A 

MSE  « 


where, 


"  8|. 


(*-«0 


p2  2 

-  R  oA 


(5-35) 


(«m  +  *,)  +  (4"ra)/b  +  (4Hrad)1 


(5-36) 


*  In  this  section  we  presume  that  Lrfts  L  in  that  the  test  is  conducted  with  flights  no 
oloser  than  8- to  12-hours  apart. 
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It  is  of  interest  to  digress  a  moment  to  discuss  the  random  error.  It  is  easy  to  show  that 
aircraft  with  mean  flight  paths  that  are  cartesian  linear,  but  not  overflying  the  radar,  can 
have  equivalent  polynomials  in  radar  coordinate  space  that  are  quite  high,  especially  at 
dose  range.  Reliance  upon  simple  flight  path  models,  without  a  precision  standard  refer¬ 
ence  radar,  can  thus  be  shown  to  lead  to  inordinately  large  model  errors  that  overwhelm 
the  radar  error.  Here  we  are  adequately  correcting  the  altimeter  in  a  few  range  intervals 
and  any  deviations  from  constant-height  flight  will  show  up  identically  In  HR  and  Hc.  That 
is  why  we  can  estimate  the  total  mean-squared  error  here  with  such  confidence,  hi  effect, 
we  have  created  a  precision  reference  standard. 

Introducing  Equation  (5-36)  into  Equation  (5-37)  and  noting  that  all  the  errors  are 
N(0,  o)  and  uncorrelated, 

E(MSE)  =  C2  +  R2  (g  ♦  <r2)  ♦  ^  +  <5-! 

which  contains  all  the  elements  we  wish  to  threshold,  i.  e. ,  the  point-to-point  random  error 
component  o  ,  the  leg-to-leg  random  bias  component  ,  the  altimeter  error  after  radio¬ 
sonde  correction  and  any  fixed  bias  term  which  may  or  may  not  be  present.  What  is  the 
error  distribution  of  MSE?  The  distribution  of  the  height  residual  Hin  is  N  (C,  <rg) 
where 


afl 


2  -  R2 


+  °RA  +  ffRAD  * 


(5-38) 


Since 


MSE 


LN  iC  S  “in  ■  °A  • 


(6-39) 


f  n 


N 


the  distribution  of  ^  n/aj|  *s  noncentral  chi-square  withN-degrees  of  freedom. 

n 

Thus, 


MSE 


o~ 


T&-  £  <N.»  -  S2<rA 


(6-40) 
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where 


X  4 


N 

L 


1  ("a  )l/4  =NC  V4- 


(5-41) 


The  non  central  chi-square  random  variable  which  arises  from  the  sum  of  the  squares  of 
non  zero-mean  Gaussian  random  variables  with  the  same  variance  is  a  difficult  distribution 
to  work  with.  A  common  practice  is  to  approximate  it  with  a  chi-square  variable  with 
appropriately  adjusted  moments  (see  reference  6,  page  942).  Thus 


where 


and 


#2  (N,  X  )  -  xl  *  d+b). 


b  =  A/(N+A), 
a  =  N+A, 

Ii  B  h  - 

v  1+b 


-.(NtA)j, 

N  +  2  A) 


(5-42) 

(5-43) 

(5-44) 

(5-45) 


Note  that, 


,2  <*> 

CN 


f-1  -t/2 

t  e 

N/2 

2  r  (N/2) 


(5-46) 


with  r  (x)  the  gamma  function, 


r  (X)  -  j 
0 


t  X-1  e_tdt . 


(6-47) 


Now  the  sum  of  L  chi-square  random  variables  with  v  *  degrees  of  freedom  is  still  chi- 
square  but  with  L  v  *  degrees  of  freedom.  Thus, 


MSE  » 


(14b)  2  2 

-  it  cr. 
LN  A 


(6-48) 


Another  way  to  look  at  this  equation  is  to  say  that  the  random  variable 


'  LN 

(1+b)  <r| 


^MSE  +  R2 


(5-49) 


is  approximately  chi-square  distributed  with  L  v*  degrees  of  freedom.  What  then  is  the 
distribution  of  MSE  ?  Equation  (5-48)  is  of  the  form 


y  -  V  -  k2 

where  x  is  chi-square.  Then 


fY(y)dy  =  £j^(x)dx 


and 


fy(y)  ■  ** 


or 


MSE 


LN 

(y)  - 


LN(y+R2  a2  ) 

L  (1+b)  o'! 


Lv* 


-  1 


2 


LN 


L  (1+b)  og 


(1+b)  <r|  2  Ll,*/2  T(LvV2) 


(5-60) 


Before  making  risk  computations,  it  seems  appropriate  to  check  on  the  adequacy  of 

the  approximate  formula  for  the  noncentral  chi-square  distribution,  i.  e. ,  Equations  (5-42) 

through  (5-45).  A  mathematical  experiment  was  performed  whereby  the  squares  of  ten 

samples  from  a  Gaussian  random  number  generator  (c»  10)  were  summed  many  times  (i.  e. 

10000)  and  then  histogrammed.  In  one  case  the  Gaussian  distribution  was  zero-mean,  i.  e. , 

N  (0,  10),  as  a  control.  In  the  second  case,  the  distribution  was  N  (10,  10).  Here  we  are 

using  v  m  io,  the  lowest  number  of  degrees  -of-freedom  of  Interest  to  us  (about  one-range 

interval  worth).  The  empirical  cumulative  distributions  are  plotted  in  Figure  5-16  as  tri- 

2  2  ° 
angles  and  circles  for  X  and  *  respectively.  The  theoretical  distributions  are  plotted 
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as  the  solid  line  running  nearly  perfectly  through  the  data.  In  the  #2  case,  the  approxi¬ 
mation  is  used  for  the  theoretical  model.  For  the  upper  tail,  I  have  also  shown  the  90% 
upper  and  lower  confidence  bounds  for  10000  Monte-Carlos.  Note  how  well  the  theoretical 
models  fit  the  data  from  this  point  of  view. 


6. 6.  3  PRODUCERS  RISK 


A  A 

Unlike  the  distribution  of  c  -  c  ,  which  was  zero-mean  Gaussian  with  variance  crS  . 

A  C 

the  distribution  of  MSE  is  much  more  complex,  being  related  to  the  chi-square  density 
function  in  Equation  (5-50).  But  the  formulation  of  risks  is  similar.  The  producer's  risk 
of  falsely  failing  due  to  inadequate  sampling  is 


RMSE  -  Pr  [mSE  >  MSEp/p/MSE  *  MSEp/Fi  |o 


where 


is  the  true  mean  square  error, 


is  the  sample  estimate  of  MSE, 


MSEp^p  is  the  pass/fail  threshold,  and 


is  the  magnitude  of  the  fixed  bias  error. 


Thus, 


■/ 


f  a  <y)dy 
MSE  * 


where  f*V  (y)  is  given  by  Equation  (5-50)  with  |  c  |  *  0.  Thus, 


LN(y+*2°A 


exP  "  T 


1  LNfy* 


MSE  (yi  1  c  !=0)=  ~ 


r  (LN/2) 


(5-51) 


(5-52) 


Transforming  variables,  i.  e. , 


(5-53) 


RP, 

nMSE 


/  /  x  fx2  (t)dt 

LN  (mSEp/f  a2  a*  )  LN 


(5-M) 


where  f^2  (t)  is  just  the  chi-square  distribution  given  in  Equation  (5-46). 

5.6.4  PASS/FAIL  THRESHOLDS 

p 

From  Equation  (5-54)  we  see  a  one-to-one  relationship  between  MSEp^F  and  R 

o  -  . 

with  parameters  LN,  a  g ,  and  R  .  The  pass/fail  threshold  for  ^MSEp^F  is  plotted  in 

Figure  5-17  for  the  same  simulated  test  parameters  as  for  the  first  test  type,  i.e. ,  L=10, 

N  =  7,  and  the  error  budget  previously  described.  The  top  two  curves  are  shown  for  risks 
of  0. 15%  and  10%  respectively.  The  bottom  curve  is  a  plot  of  erg  versus  range  showing  how 
linear  the  pass/fail  thresholds  are  in  this  parameter  even  with  the  complexity  of  Equation 
(5-54).  It  is  important  to  reflect  here  that  these  computations  imply  radiosonde  correction. 
Even  though  the  root-MSE  pass/fail  threshold  can  be  many  thousands  of  feet  for  ranges 
beyond  120  nmi,  it  appears  to  be  relatively  insensitive  to  risk  level. 


5. 6.  5  BUYER'S  RISK 

Similarly  to  Equation  (5-27), 


RMSE  "  pr  35  ^Ep/F  / mE  >  mEP/F  ’’  iC  lCB  (6'*56) 

where  |Cfi  |  is  the  maximum  fixed  bias  error  the  buyer  must  tolerate  due  to  sampling  in¬ 
adequacy,  Thus, 
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Pass/Fall  Threshold  for  the  Total  Mean  Square 
Height  Error  Test 


A  A 

MSE 


/ 


MSE 


P/F 


*MSE  ^  y*  1CB  1  )  ^y* 


(5-66) 


•  5  2  2 
-R  <rA 


We  will  treat  all  negative  estimates  of  MSE  (possible  because  of  the  correction  factor 
—  o  2 

-  R  <rA  )  as  zero.  This  is  equivalent  to  considering  only  positive  values  of  the  distribution 
with  an  impulse  at  zero  equivalent  to  the  negative  area.  Thus. 


R®~ 

MSE 


MSE-d  /it  ( 


*me  (y-  l°B  !)  *  (t-BV 


B 


-R  ai 


where  the  second  term  is  the  power  of  the  Impulse. 
Making  the  variable  transformation 


t  A  LN 


(1-Ha)  o| 


(y+*2°l ) 


we  have, 


(MSEp/r*E2<ri  )  (1^  2 

J  fx?u*(t)dt  +  J  £v2  (t)dtf 


LNR2oA 


(5-58) 


.B 

l  A 

MSE 


with 


LNR2a2 
(1+b)  o| 


X  =  NC|/<r|  . 
a  -  N+X  , 
b  -  X/(N+A)  , 


V  * 
Lv* 


(6-59) 

(5-60) 

(5-61) 

(5-62) 
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v*  «  a/(l+b)  -  (N+A)  /(N+2X), 


(6- 

Note  that  the  second  term  has  the  same  value  whether  we  treat  negative  values  of  MSE  as 
zero  or  use  them  as  they  are.  Thus  we  have  that, 


XB 

R^E  "  /  V  <4>dt 

J  *Lj>* 

with 


LN 

(1+b)  <r| 


(y+x2°Z)  ■ 


(5- 


(5- 


Solving  Equations  (5-59)  thru  (5-65)  for  Xg,  b  and  Lv*  in  terms  of  (Cfi  |  ,  we  get 


Lv* 


2Cb+4  4 

(4+2Cb)4 


and 

b  -  c! /(4  +  cl)  ■ 


(5- 


(5- 


(5- 


B 

These  equations  provide  a  unique  one-to-one  relationship  between  j  |  and  Rj^e  with 
parameters  LN,  ,  and  MSEp/-. 


5. 6. 6  BUYER’S  LIMIT  FOR  C„ 

Xj 

Equations  (5-64)  thru  (5-67)  are  plotted  in  Figure  5-18  for  a  producer’s  risk  of  10% 
and  a  buyer's  risk  of  10%  overall.  The  topmost  curve  is  the  standard  deviation  of  the  height 
residual,  ag  .  The  other  three  curves  are  the  buyer’s  limits  on  C  for  (rp«  10%,  rB«40%), 

rB  “  96%),  and  (rP«  io%,  rB«80%)  respectively  corresponding  to  (M-l,  N  «5), 
(M  =  0,  Nq  «  70),  and  (M  =  10,  Ns«  70)  as  before.  Note  that  {Cfi  |  is  larger  in  the  MSE  test 
than  in  the  fixed  height  bias  test.  This  is  a  consequence  of  summing  over  squared  residuals 
instead  of  the  residuals  themselves  as  in  the  fixed  height  bias  test.  The  comparison  between 
the  bottom  two  curves  of  Figure  5-18  is  interesting.  Here  we  see  the  material  benefit  to  the 
buyer  of  permitting  failures.  At  120  nm*  for  instance  the  M=0  curve  yields  a  buyer’s  limit 
of  620  feet.  Permitting  10  out  of  70  failures  and  allowing  the  producer  to  still  pass  the  test 
reduces  his  limit  on  the  fixed  height  bias  of  380  feet,  almost  a  factor  of  2. 


5. 6.  7  BUYER’S  LIMIT  ON  MSE 

From  Equation  (6-48)  we  have  a  general  relationship  between  estimated  MSE  and 
the  chi-square  random  variable.  Thus, 
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Figure  5-18.  Warimmn  Tolerable  Bayer's  Limit  on  Fixed  Hei 
Bias  Error  as  Part  of  a  Total  Mean  Square  Heij 
Error  Test 


This  is  plotted  in  Figure  5-19  for  the  same  pass/fail  criteria  of  Figure  5-18.  Onoe  again, 
the  curve  for  eg  is  plotted  (as  a  dashed  line  here)  for  comparison.  Again,  radiosonde 
data  is  presumed  available  to  correct  the  Mode-C  height.  In  this  instance,  unlike  for  the 
buyer's  limit  on  the  fixed  component,  1.  e. ,  j  j ,  there  is  only  marginal  benefit  of  per¬ 
mitting  failures  insofar  as  the  buyer's  limit  on  MSE  is  concerned.  Nevertheless  the  buyer 
gains  overall  by  forcing  the  producer's  risk  per  subtest  up  from  0. 15%  to  10%  and  by  re¬ 
ducing  his  risk  per  subtest  down  from  96. 5%  to  C0%.  His  limit  on  the  fixed  component  of 
height  bias  error  is  of  course  materially  reduced  as  well. 
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